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^L ■ The p-adic formulation of replica symmetry breaking is presented. In this approach 

\^ . ultrametricity is a natural consequence of the basic properties of the p-adic numbers. 

^^ I Many properties can be simply derived in this approach and p-adic Fourier transform 

Qv ' seems to be an promising tool. 
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1 Introduction 

In the replica approach to disordered systems one usually introduces a matrix Qa,h which is 
the stationary point of a free energy -F[(5]; the matrix is a zero by zero matrix, with zero 
elements on the diagonal [1]. Such a matrix is constructed as the n — *> limit of a normal 
matrix with n components. 

In the mean field approach one looks for stable (or marginally stable) saddle points of 
the free energy. When the replica symmetry is spontaneously broken, as it happens in spin 
glasses, one assumes that the saddle point is given by a matrix Q constructed in a hierarchical 
way, which corresponds to breaking the replica symmetry group (the permutation group of 
n elements) in a peculiar way [2, 3]. The aim of this note is to expose some hidden algebraic 
properties of this matrix and to show that the whole construction may be simply done using 
p-adic numbers. 

In this approach the ultrametric properties of the matrix Q [4, 5] arise naturally from 
the ultrametric properties of the p-adic numbers. Although we do not obtain new results in 
this way, we hope that this reformulation may be a useful starting point for simplify some 
of the long computations involved in the evaluation of the corrections to the saddle point 
approximation. 

In section 2 will be present the basic properties of the p-adic construction and show its 
equivalence to the usual hierarchical construction. In the next section the limit n — ;► is 
performed in a simple way. In section 4 we present an alternative and more interesting 
procedure for doing the limit n ^ 0, where we connect this approach to standard p-adic 
analysis. In the next section we show the advantages of using the p-adic Fourier transform. 
Finally there are two appendices, the first dedicated to the foundations of p-adic analysis, 
the second to the basic properties of the Fourier transform[6, 7]. Both appendices can be 
skipped by readers experts on p-adic analysis. 

2 The p-adic construction of the matrix Q 

We start the construction of the matrix Q by considering a number p (which for simplicity 
we suppose to be a prime number) and by assuming that n = p^ for some value of L. We are 
going to construct the matrix Q for integer L and p in a specific way which we will discuss 
later. The limit n — > will be done at the end. The matrix Q enters in the evaluation of 
the free energy in spin glasses and related models in the saddle point approximation. Here 
we do not address the point of the evaluation of the free energy and we only consider the 
construction of the matrix Q. 

Eventually n = p^ must go to zero and we can follow two options in order to realise this 
goal: 

a) We take a value of p greater than one and we send L to minus infinity. Eventually we 
may do an analytic continuation in p to non integer values of p. 

b) We first do an analytic continuation in p to non integer values of p. We take a value 
of p less than one and we send L to plus infinity. 

In both cases one obtain the limit n ^ 0. The two constructions are roughly equivalent. 
It seems that the second one is more simple to work with, however for pedagogical reasons 



we will start by presenting the first one in section 3, while the second one will presented in 
section 4. 

The first steps are common to both strategies. The construction of the matrix Q for 
integer p and L can be done as follows. We assume that the matrix Qa.t is of the form 

Qa,b = Q{a-b). (1) 

where Q{k) = Q{—k) (symmetric matrix) and Q{k + n) = Q{k). This choice restricts very 
much the form of the matrix and shows an explicit symmetry of this parametrization (i.e. 
translational invariance in internal space). The condition Q{0) = implies that the elements 
on the diagonal are equal to zero. 

The second step consists in assuming that the function Q{k) is a function of the p-adic 
norm \k\p. The appendices provide a brief introduction to peadic analysis. 

In other words we suppose that 

Q{k) = q{\k\,). (2) 

This corresponds to setting 

Qa,b = qi = q{p~'), if \a-b\p = p~\ (3) 

Before performing the limit n — *> it is convenient to compare our apporach with the 
standard hierarchical construction. 

In the usual case [2] one introduces a sequence oi K + 2 numbers rrii, with mo = 1 and 
rrix+i = n, such that mj_i divides rrii for i = 1,K + 1. One sets for a ^ b: 

Qa,b = Qi if lio-l^i) 7^ Hb/rrii) and /(a/mj+i) = /(6/mj+i), (4) 

where the function I{z) is the integer part of z, i.e. the largest integer less or equal to z. 
Let us consider the special case where the rrii are given by 

mi=p\ (5) 

We want to show that the matrix Q obtained in this way coincides, after a permutation with 
the matrix Q constructed before with K + 1 = L. The proof is rather simple. We associate 
to the index a the L digits of a — 1 in base p: 

a = l+ J2 «^P' (6) 

i=0,K 

These digits form a L dimensional vector with components in the range — {p — 1). The 
hierarchical construction corresponds to set Qa,b = Qi if dj = bj for all j > i and aj_i ^ 6j_i. 
We now associate to an index a its transpose (ot), which is obtained by writing its digits 
in the inverse order: 

ar = l+ Y: a,p(^-). (7) 

i=0,K 

The previous condition becomes that the K — i less significative digits of a^ and b^ do 
coincide, and the {K — i + ly^ digit differs. This last condition may be restated by saying 



that a — 6 is a multiple of p(^~*) but not of p(^^*+^), i.e. \a — b\p = p~(^"*). Apart from a 
reshuffling of the indices, i.e. a permutation, the usual construction is equivalent the p-adic 
construction presented before. 

Generally speaking it is possible to prove that independently from the condition in eq. 
(5), after a similar reshuffling of the indices, the hierarchical matrix Q (defined in eq. (4)) 
can always be written under the form of eq. (1). It is likely that many of the unexpected 
properties of the hierarchical construction arise from the possibility of choosing an ordering of 
the indices in such a way that the hierarchical matrix is invariant of under the transformation 
a ^ a + 1. This invariance implies that the elements of one line are the permutation of the 
elements of another line, however the converse is not true. 

3 The n ^ limit 

We can now perform the n ^ limit. We will firstly follow the strategy a). 

This limit can be reached by sending L to — oo. The continuation of the usual formulae 
from positive to negative L can be done if we introduce the quantities qi for non positive i. 

For example let us consider the sum of the elements of a line of the matrix. In order to 
perform the n — > limit we slightly change the notation of the previous section and we set 

Qa,b = Qi, if \a - b\p = p^-\ (8) 

We easily get that the sum is given 

T.Q<^,b= E(p-iK"'^- (9) 

b i=l,L 

Indeed the number of integers k such that |A;|p < p~^ (i.e. the volume of the p-adic sphere) 
is given by p^~^ and therefore the number of integers k such that |A;|p = p~^ (i.e. the volume 
of the p-adic shell) is given by {p — l)p^~-'~^. 
We are free to write the last equation it as 

E ip-'^)p''%- E ip-^)p'~'<l^, (10) 

i=—OD,L i=— 00,0 

by introducing the extra parameters q^ for i < 1, which are irrelevant for positive L. 

The analytic continuation to negative L can be now trivially done. In the limit L ^ — oo, 
the first term disappears and we get 

EQa,6^ E ip'''-p%- (11) 

b i=—OD,0 

A similar procedure can be followed in order to compute other functions of the matrix 
Q. By comparing the previous equation with the usual ones, we see that we obtain the 
hierarchical formulation where a function q{x) is introduced, with the extra constraint that 
q{x) is piecewise constant with discontinuities at a; = p~*. The usual formulation, where q{x) 
is a continuous function can be obtained by analytic continuation in p up to the point p = 1^. 

By performing the explicit computations similar results are obtained for the other quan- 
tities and it is possible to show that the usual approach is recovered. 



4 The upsidedown world 

In the other possible approach to the n — > hmit (b), we firstly do an analytic continuation 
in p to values less than one and only later we send L -^ +cxd in such a way that p^ -^ 0. At 
a later stage we are free to send p ^ 1~ in order to reach the continuous limit. In this way 
we get formulae quite similar to the previous one, with the advantage that only the qi with 
positive i are needed. 

In this approach one obtains that the function q{x) is given by 

q{p') = q^, (12) 

where the index i ranges from to +oo in such a way that when p ^ 1~ , x = p^ spans the 
interval 0-1. The formulae one obtains in this approach for p < 1 coincide with the formulae 
obtained with the formalism of the previous section (with the substitution of p with p~^). 

The advantage of this procedure is that we obtain formulae that are very similar to those 
used in the p-adic integral and that are well known to mathematician. The strategy to prove 
these formulae is quite similar and therefore one can use some of the well known results in 
this field. 

In the region where p < 1 it may be convenient to introduce the notation: 



\k\ 



\k\ = ^. (13) 

In this way \k\ belongs to the interval — 1, with the exception |0| = oo. In the limit where 
p —* 1~ , \k\ spans the interval 0-1. Equation 12 can thus written as 

Q{a-b)=qi\a-b\). (14) 

Let us apply this strategy to the computation of the sum of the elements of a line of the 
matrix. We find that 

hm J2 Q(«-&)= E(P-iy"'?- (15) 

L-*oo — ^ . ~^ 

a=l,p^ j=l,oo 

For p < 1 the previous equation can be written as 

j2Q{a-b)= Y. a-p)p'-'q^, (16) 

a=l,n i=l,OD 

while for p > 1 the r.h.s. becomes proportional to the p-adic integral which is denoted as 

daQ{a). (17) 



'p 
With some abuse of notation we denote for p < 1 



lim y Q(a-b)= f daQ(a), (It 

'^^0 „ 1 „ Jp 



a=l,n 



where the sign ' over the integral / denotes that the value zero is excluded from the integra- 
tion range. We must note that the measure of the integral is normalised to —1. In a similar 
way we can use the notation 

lim y~] F{\a — b\,\b — c\,\c — a\) = db daF{\a — b\,\b — c\,\c — a\). (19) 

^^^ u 1 L h 

For p > 1 we obtain the usual p-adic integral (apart from a normalisation factor). The results 
for p < 1 can be obtained using the same steps as in appendix II. 

We can do the computation in the interesting case where the sum is restricted to all 
different indices. We have to compute 

lim V" F(|a — 6|, |6 — c|, |c — a|) = 

lim V 'F{\a-bl\b-cl\c-a\)= (20) 

b,c=l,p^ 

f dbdcF{\a-b\,\b-c\,\c-a\)= (21) 

Jp b,c=l,p^;ay^b,ay^c,by^c 

/ da dbF{\a — b\,\b — c\,\c — a\), 
Jp 

where we denote by J2 the sum restricted to the case of all different indices. 
In the same way we could define 

1 /■' 

(-l)^lim- y F\a]= daida2...daM,F[a] (22) 



ai,a2,...aM 



where p is less than 1 and the sum is done on all different indices. The factor (—1)^ has 
the effect of giving a positive result for the integral. Generally speaking in this way the 
evaluation of sums can be reduced to the computation of quantities that are very similar to 
the corresponding p-adic integral. 

For example let us use this strategy to compute 

f dadbF{\a-b\,\b-c\,\c-a\). (23) 

Jp 

The application of the previous formulae tells us that the integral is given by 
E F{\a-b\,\b-cl\c-a\) = EM^, K l)F{p',p\p^) = 

b,c,bj^a,Cy^a,by^c i,k,l 

E ip'^' - P'){p'^' - p')[F{p',p',p') + F{p',p\p'') + F{p\p',p')] + 

i,k\i<k 

E(P^+' -pW^' - 2f)F{p\p\p^). (24) 

i 

The proof can be obtained using the same strategy as in the appendix I for computing 
the measure of three intersecting p-adic shells. 



Finally in the continuum limit where p goes to 1 one get the formula 

^ F(\a-b\,\b-c\,\c-a\) 

b,c,by^a,Cy^a,by^c 



dxdy6{x — y)[F{x,x,y) + F{x,y,x) + F{y,x,x)] + dx x F{x,x,x). (25) 

The same formula could be simply written as 

y^ F{\a — b\,\b — c\,\c — a\) = dxdydxfj,{x,y,x)F{x,y, z) 

b,c,byta,Cy^a,by^c 

fi{x, y, z) = 9{x — y)6{x — z) + 9{x — z)6{x — y) + 6{y — x)5{y — z) + x5{x — y)S{x — z) (26) 

It is important to note that the ultrametricity inequality works at reverse in the region 
p < 1 and consequently also in the limit p ^ 1^. This is in agreement with the fact that 
1 — x, not X, has the physical meaning of distance. 

After some work one can find simple rules for generic sums of the type 

— ^ 'F{a,b,c,d)= / F{a,b,c,d) (27) 

^ a,b,c,d •'v 

where F depends only on the p-adic distance and all indices are different [9]. In the case 
where the function is symmetric one finds that in the limit p ^ 1^ 

Ip F{a, b, c, d) = 

Ix<y<zdx dy dzF\ia--b\=\a-c\=\a-d\=x,\b~c\=\b-d\=y,\c~d\=z+ H permutations + 

L<y;x<zdx dy dzF\ia-b\=z,\b^c\=\b-d\=\a~c\=\a-d\=x,\c-d\=y + 2 permutatioHs 

+ lx<ydx dy X F\\a-b\=y,\b-c\=\b-d\=\a-c\=\a-d\=x=\c-d\=x + 5 permutations 

+ L<yd^ dy y F\\a-b\=\b-c\=\a-c\=y,\b-d\=\a-d\=x=\c-d\=x + 3 permutations 

+2\ J dx X F{la-b\ = \a-c\ = \a-d\ = \b-c\ = \b-d\ = \c-d\=x), (28) 

where the formula for the intersection of 4-p-adic sphere of the same radius has been crucial 
to obtain the last term. 

If we apply the same strategy to more complicated sums we can find the formula of ref. 
[8], where the result is written as sum over all possible trees, with a specific integral associated 
to a given tree. 

5 Using the p-adic Fourier Transform 

An interesting application of this approach can be done to the formula for the product of 
two matrices A and B: 

j 
If the matrices have the form discussed here one finds that the previous formula can be 
written as a convolution 

C{t - k) = Y.A{i-j)B{j - k) = Y^a{\i-3\)b{\j - k\). (30) 



Finally one finds using the previous formulae that 

c{\i\) = f dka{\i-k\)b{\k\), (31) 



which using the rules of p-adic integral after performing the limit p ^ 1 can be written as: 
c{x) = / dy{a{y)b{x) + a{x)b{y)) + / dya{y)h{y) + x a{x)h{x). (32) 

J X J X 

Convolutions can be strongly simplified in Fourier space. In principle we can just do ordinary 
Fourier transform, where the momentum q is in the interval (— tt , tt), however it is convenient 
to take into account the p-adic nature of the functions we consider. 

We can start from the analysis leading to formula (64) of appendix (II). Generalising the 
computations to the case where p < 1 and performing the continuum limit, one finds that 
the p-adic Fourier transform^ of a function A{k) = a{\k\) is a function A[M] = a[y], defined 
for y in the interval 0-1, where y = \M\^^. One must be careful in removing the factor ^ in 
the normalization of the Fourier transform, which would be harmful here. 

Eq. (64) may be transformed to 

A[M]||M|,=p. = E (^"'^"^(1 - l/pW-k) - p-'^'-aL-, + A(0). (33) 

fc=0,i 

Let us define (for p < 1) the function 

a[y]=A[M]\\M\,=r,Ly (34) 

The function a[y] is defined only for y of the form p~^ for integer k. With this definition one 
finally obtains that 

a[y] = J2 /(I - ^/pHp~') -P^a{p~^) + A(0). (35) 

fc=0,i 

The final formulae in the continuum limits are 

a[y\ = a(oo) — / dxa{x) — ya{y), 

a\oo\ = a(oo) — / dxa{x), (36) 

Jo 

where we use the apparently strange notation a(oo) = A{0) and a[oo] = A[0]. 
The inverse Fourier transform formulae simply read 



a{x) = a\oo] — a[0] + / 

Jo 



dy da[y] 



y dy 

a(oo) = a[l]. (37) 



^We use the same notation as the appendix and we use the square parenthesis, [•], to denote the Fourier 
transform 



These relations becomes simpler in differential form. For example one gets: 

X a'{x) = a'[x]. (38) 

This differential relation is equivalent to the integral relation in eq. (36, 37) for obtain the 
Fourier transform if they are complemented by value of the Fourier transform in a given 
point. A possible choice is 

a[l] = a(oo) — a(l). 

With some work one can verify that the multiplication of two matrices becomes the simple 
multiplication of their Fourier transform [10]: 

c[y] = a[y] b[y]. (39) 

Indeed diffentiating eq. (32) one finds 

c\y) = a'{x){[ dyh{y) + xh' {x)) + h' {x){ ( dya{y) + xa {x)) (40) 



The function a[y] was already introduced in ref. [10] in order to solve the inversion prob- 
lem, although its p-adic nature was not recognised. Also the usual procedure of simplifying 
the saddle point equations by differentiating them correspond to consider the p-adic Fourier 
transform. 

We could also consider the problem of computing the inverse of a matrix Q, i.e. of finding 
matrix R, such that 

Y. Qa,cRc,b = Sa,b. (41) 

c 

It is not a surprise that we find the Fourier transform of the matrix R is simply the inverse 
of the Fourier transform of the matrix Q: 

q[y] = l/r[y]. (42) 

An extremely important problem consists in the computation of the inverse of the Hessian 
coming from the fluctuation around the saddle point. Here one has to solve the equation 

/ , Mafi-,ejGeJ;c,d = ^a,b;c4- (43) 

e,/ 

This inversion is not a simple job and rather complex computations have been done [11, 12]. 
However the final formulae are remarkable simple. Although we are not able at the present 
moment to derive these formulae in the framework of the p-adic formalism it may be useful 
to show that they have a very simple interpretation in terms of p-adic Fourier transform [7]. 
We will consider here only the so called replicon sector for which the results are simpler 
than in the other sectors. We restrict our analysis to the region where \a — h\=z>\a — c\ = 
Xi, \h — d\ = X2 > z. In this region ultrametricity implies that z = \a — d\ = \b — c\ = \c — d\. 
Both M and G are functions of a;i, X2, z only and we write them as M^{xi,X2) and G^{xi,X2)- 
In the same way we denote by Gf^{xi,X2) the replicon contribution to G, where the precise 
definition of the replicon can be found in the original papers [11, 12]. 

9 



Following [11] we can thus introduce the Fourier transform with respect to Xi and X2, 
which is given by 

nf{xi, 1/2] = ^^{xi, 00) - / dx2 'rrf{xi,X2) - y2'm^{x 1,1/2), 

m^[yi,y2]=m^{oo,X2]- dxi m''{xi,y2\ - yim^{xi,y2]. (44) 

■Jyi 

One finally finds that the final formula for the replicon propagator [12] may be obtained 
with slightly modified inverse Fourier Transform: 

52 92 1 

9a;i9a;2 ' 8x18x2 m^[xi,X2\'' 

9r{z,X2) =g%{xi,z) =0. 

Equivalently we could write the last equation as 

rz f^ I 8^ 1 
gR{xi,X2)= dxi dx2 ^-^ -, 7. (46) 

Jxi Jx2 X1X2 8X18X2 m^[Xi,X2\ 

The differential relations in the inverse Fourier transform are preserved, only the second 
condition which fixes the value of the inverse Fourier transform in one point is modified. 
With these modifications the replicon sector of the inverse is just the numerical inverse of 
the matrix M in Fourier space. 

The precise reason for the appearance of these simple formulae with a strong p-adic 
fiavour is not completely clear at the present moment. They show the usefulness of the p- 
adic formalism. It would be also extremely interesting to study if the same formalism could 
be applied to the off equilibrium dynamics of the kind studied in ref [13]. 

After completion of this work we received a paaper by V. A. Avetisov, A. H. Bikulov, S. 
V. Kozyrev [14] where some similar results are derived. 

Appendix I: p-adic numbers 

Let us consider a prime number p. Any integer k can be written in an unique way as 

00 

k = p'J2(^iP' (47) 

1=0 

with i>0, 0<a;<p — 1 and ag 7^ 0. The p-adic norm of such an integer k (i.e. \k\p) is 
defined as 

\k\,=p-\ (48) 

The p-adic norm of is defined to be equal to zero. The value of the p-adic norm tells us 
the number of consecutive zeros at the end of a number, when it is written in base p. For a 
rational number r = a/b, the p-adic norm is defined as \r\p = |a|p/|6|p. 

10 



The properties of the p-adic norm are well studied by mathematicians, one of the most 
famous property being ultrametricity, which states that 

|o — b\p < max(|a — c\p, \c — b\p) (49) 

for any choice of c. This property, which generalises the statement the sum of two even 
number is even, can be proved as follows. 

Using the translational invariance of the metric, we first write the ultrametric inequality 
in an equivalent way as 

\a + b\p < max{\a\p,\b\p). (50) 

If a is a multiple of p* (and not of p^^^), and 6 is a multiple of p'', with k > i, it is evident 
that a + 6 is a multiple of p^. Therefore 

|a + b\p < p~^ = \a\p = max(|a|p, \b\p). (51) 

We stress that we have used in a crucial way the fact that p > 1 for a true prime. (In this 
paper we make an analytic continuation to p < 1. In that case, the inequality sign would be 
reversed.) A direct consequence of this inequality is that any triangle is either equilateral or 
isosceles with the two largest sides equal. It follows that any point a inside the p-adic disk 
centered at o and of radius r, i.e. such that \a — o\p < r, is also a center of the disk; i.e. if 
1^ — o\p < T, then also \a — b\p < r. 

The whole p-adic field may be constructed starting from the p-adic rationals by consid- 
ering the closure of the rationals with respect to the p-adic norm in the same way that the 
real numbers (of the interval — 1) are constructed as the closure of the rationals (of the 
interval — 1) with respect to the usual Euclidean norm. 

Closing the rational field with respect to the previously defined norm one obtains the 
p-adic field. Continuity of a p-adic function can be defined as usual. For example a function 
/ is continuous at the point k if 

jim/(A;„) = /(A;), (52) 

for any sequence of A;„ which converges to k in p-adic sense (i.e. \kn — k\p -^ 0). The extension 
of a function from integers to p-adic numbers is called p-adic interpolation. Here we do not 
need to discuss this point any more. 

/^From our point of view a more interesting construction is the integral over the p-adic 
integers which can be defined in an elementary way as 

1 "' f 

lim —yF(a) = / daF(a). (53) 

L^oo pL ^^ Jp 

There are many well known properties of the p-adic integral. Here we report some of 
them, leaving the proof to the reader (the lazy reader can found them in any book on p-adic 
calculus) . 

a) The measure of the p-adic sphere of radius p^^ centred around an arbitrary point a 
(i.e. the measure of all points such that \a — b\p < p~*) is given by p~\ As far as the p-adic 
distance among integers cannot be larger than 1, the unit sphere coincides with the whole 
space and has measure 1. 

11 



b) The measure of the p-adic shell of radius p~* centred around an arbitrary point a (i.e. 
the measure of all points such that \a — h\p = p"*) is given by p~* — p^^^^ = (1 — p^^)p^^. 

c) The measure of the intersection among two p-adic shells has rather interesting prop- 
erties. Let us consider the intersection of a shell of radius p~* centred around the point a 
with a shell of radius p^'' centred around the point b. The measure depends on the distance 
among the points a and b, which we assume to be equal to p~^ . Ultrametricity tell us that 
the measure is zero unless two among the distances coincide and the two equal distances are 
the largest. After some reflection one finds that only three cases have to be considered. 

• We first consider the case i = k < j. Here the ultrametricity inequality implies that 
the two shells coincide and therefore the measure of the intersection is simply given by 
(1 — p~^)p~\ 

• We now consider the case i = j < k. Here the ultrametricity inequality implies that 
the second shells is fully contained in the first one and therefore the measure of the 
intersection is simply given by (1 — p~^)p~^ . 

• We finally consider the less trivial case isi = j = k. If one notice that the two spheres of 
radius p~* centred in a and b coincide and that the two spheres of radius p~*~^ centred 
in a and b have zero intersection one finds that the measure of the intersections of the 
two shells is given by (1 — 2p^^)p~~\ 

d) The generalisation of the previous arguments allows us to compute the measure of the 
intersection of many p-adic shells by using ultrametricity in a systematic way. The most 
significative result is that the intersection of M shells of radius p~*, whose centres are all at 
mutual distance p~^ is given by (1 — Mp~^)p~^. The measure becomes zero for p = M, which 
implies that you cannot find M + 1 numbers exactly at the same distance. This last results 
is a generalisation of the well known statement that you cannot find three integers (a, b, c) 
such that the three differences among them {a — b,b — c,c — a) are all odd. 

Using the previous formulae there are a few p-adic integrals that can be obtained a simple 
way. 

For example let us try to compute 

/ da dbF{\a — b\,\b — c\,\c — a|), (54) 

Jp 

where for simplicity we denote by |a| the p-adic norm of a. 

The integral is c independent and the application of the previous formulae tells us that 
the integral is given by 

j:i,j,i^jP~"P'^{l - P'yF{p-^,p-\p-^) + 

E^,J>^p-'p'''{l-p~'?F{p-\p-\p-') + (55) 

EiP"^*(l-2p-i)(l-p-i)F(p-\p-\p-^), (56) 

where k = min{i,j) 
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Appendix II: p-adic Fourier transform 

Fourier transform on the p-adic integers coincides with the usual Fourier transform. It can 
also be defined by analysing the characters of the addictive group. It is more simple to 
consider first the case where L is finite and only a finite number of points is present. 

We start by considering the case in which the function A{x) is defined only for x = 
1,- ■ ■ ,p^ (with A{0) = A{p^)). The Fourier transform is defined as 

1 ^" 
A[M] = — Yl exp(27rimM)A(m), (57) 

P m=\ 

where M is a rational number of the form jp^^ with < j < p^. As usual the Fourier space 
contains the same number of points of the original space. In this paper we will use the square 
parenthesis to denote Fourier transform. 

Let us consider the problem of computing the Fourier transform of a function which 
depends only on the p-adic norm,i.e. A{k) = a{\k\p). We are thus interested in computing 

1 P^ ^-^ 1 

A[M] = — Yl exp(27ri mM) a{\m\p) = J2 a{p-'')Sk[M] + — A(0), (58) 

P m=\ fc=0 P 

where S\^ [M] is the Fourier transform of the p-adic shell of radius p 



-k. 



P^ 



Sk[M] = — J2 V|,P-* exp(27rz mM). (59) 

P m=l 

In order to compute the Fourier transform of the shell, it may be simpler to firstly compute 
the Fourier transform of the p-adic sphere of radius p~^. A simple computation shows that 

Vk[M] = — y e(\m\ -p~^) exp(27ri mM) = ^ V exp(27ri mM/) = 

■r m=l " m=l 

exp[27rtlp'' ) — I 

where we used the definition M = lp~^ . 

It follows that Vfc[M] = unless / = np^~'', n = 0, 1, ■ ■ ■ ,p — 1 and the last passage is no 
more valid, because both numerator and denominator are equal to zero in the final result. We 
notice that the possible values of |M| are / for non negative j. Consequently we find that 
\4[M] = unless k — j < We also remark that, as consequence of translational invariance, 
the Fourier transform of a function of the p-adic norm, is still a function of the p-adic norm. 
We can thus define the functions a[ ] and Sk[] as 

a[p'^] = A[M] \\M\=pJ 
Skip-'] = Sk[M] ||M|=p. (61) 

The reader should notice that the functions a[ ] and Sk[ ] are defined in such a way that it 
argument is the range — 1. It follows that 
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Vk [p -'] = p '^ for k — j < 

Vk[p~^] = for k-j>0. (62) 

As a consequence we find that the Fourier transform of a spherical shell is given by 

Skip'^] = P~^(l - 1/p) for k- j <0, 

Sk[p'^] = -p-^ for A; - j = 0, (63) 

Sk[p~''] = for A; — j > 0. 
The final expression for the Fourier transform thus becomes 

a[p~^] = ^ p-\l - l/p)a{p-'')+p'^a{p-^) + ^A(O). (64) 

k=0,j P 

It is interesting to note that the last formula the dependance on L is very simple, so that 
the limit L ^ oo can be trivially done. Moreover most of the properties of the ordinary 
Fourier transform, like the theorems concerning convolutions, are still valid. 
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